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GIBBS' STATES FOR MOSER-CALOGERO POTENTIALS 

K.L. Vaninsky* 
School of Mathematics, Institute for Advance Study, 
Princeton, NJ 08940, USA 

We present two independent approaches for computing the thermodynamics for classical 
particles interacting via the Moser— Calogcro potential. Combining the results we propose 
the form of equation of state or, what is equivalent, the asymptotics of the Jacobian 
between volume elements corresponding two symplectic structures on the phase space. 

1. Introduction 

Non-equilibrium statistical mechanics studies properties of dynamical systems de- 
scribing big or infinite systems of interacting particles. The Hamiltonian structure 
of the equations allows one to define on the phase space so-called Gibbsian mea- 
sures, invariant under the dynamics. These measures are usually constructed from 
the basic Hamiltonian(=energy) and other classical conserved quantities such as the 
number of particles, integrals of the momentum and the angular momentum. They 
constitute a finite parameter family, so-called grand-canonical ensemble. At present 
there are four known ID hierarchies with similar structure of the ensemble in the 
generic case and the same number (2) of integrable exceptions, as it is explained 
below. 

Statistical mechanics of ID particles is well studied, 0. In the generic case of 
Lenard-Jones potential the Gibbs' states form 3 parameter family. In 1975-76, 
Moser and Calogero discovered two repulsive potentials U(y) — x~ 2 , sinh~ 2 x 
which are integrable by the inverse spectral transform, an. Additional Gibbs' states 
from additional integrals of the motion different from the classical ones for such 
potentials were constructed by Chulaevskii □. In the paper of Gurevich the 
exceptional character of Moser-Calogero potentials was proved. 

Statistical mechanics of the ID semi-linear wave equation with restoring force 
was studied by McKean- Vaninsky ml The phase space of the dynamical system 
is a product of two function spaces. In this phase space, a Gibbsian measure was 
constructed and the existence of flow was proved, as was the invariance of the 
measure under the flow. 

Also, two integrable equations were considered: linear Klein-Gordon and Sinh- 
Gordon. These equations have integrals of local densities as an additional conserved 
quantities of motion. These integrals produce infinitely many additional Gibbs' 
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states. Klein and Sinh-Gordon can be considered as analogs of Moser-Calogero 
potentials. 

In such integrable equations the conserved quantities prevent any Gibbs' state 
from being ergodic in the periodic case. On the entire line, .these are destroyed on 
the support of the measure. As it was proved by McKean 0, for the sinh-Gordon 
equation, the Gibbs' state constructed from the basic Hamiltonian is ergodic. It is 
believed that all other mutually singular Gibbs' states are ergodic too. 

Another field model where it is possible to carry out this programis the de- 
focusing nonlinear Schrodinger equation with the power nonlinearity, EnI3. Again, 
there are two integrable equations, the linear and the cubic cases, with additional 
Gibbs' states. 

Recently, 0, another hierarchy of so-called modified KdV equations with similar 
properties was discovered. Now all these one-dimensional hierarchies fit nicely into 
the table. 
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The first line of the table displays the Hamiltonians of these four hierarchies. 
For particles, the potential U(-) is assumed to be Lenard- Jones type. For PDE 
hierarchies F(-) is the potential of a restoring force; it is an even function increasing 
at infinity. The next three lines are generic classical conserved quantities. The NLS- 
hierarchy and m-KdV hierarchy do not have any analog of N — D of particles. The 
next line is the symbol for the Gibbs' measure; it is a three/two parameter family. 
The last line presents integrable potentials where additional conserved quantities 
produce an additional Gibbs' states. 

All three integrable nonlinear equations are so-called AKNS systems, related 
to some selfadjoint problem for the Dirac operator. These are all the nonlinear 
equations related to such spectral problem. 

Additional structure of integrability allows one to obtain detailed information 
about the Gibbs' states. In section 2 we describe the decomposition of the canonical 
measure for the cubic Schrodinger equation previously obtained in Ej'113. In section 
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3, taking such decomposition for granted for Moser-Calogero potentials, we compute 
the thermodynamic limit of the specific free energy, the average energy per particle 
and the equation of state. In the last section 4 we present an alternative approach in 
the case of the rational potential 2x~ 2 and formulate the conjecture on the explicit 
form of the equation of state and the asymptotics of the Jacobian. 

2. Gibbs' States in Action-Angle Variables. 

For the cubic Scrodinger equation, the global action-angle variables < Ik, < 
4>k < 27r, k = — 1, 0, 1, associated with the basic symplectic structure ui were 
constructed by McKean-Vaninsky La. The next step is to express the measure 

e~ H dvo\ where H = / \ij)'\ 2 + 2|i/;| 4 

Jo 

in such coordinates. The main ingredient is the trace formula H — ^ I' k , where the 
I's are periods of some Abelian differential on the Riemann surface. Presumably 
they are the actions, but relative to some higher symplectic structure to' . On a 
formal level, 

e - H dvol = e~T, 1 ' Y[ dI # = e ~ ^ r Jac Y[dl'dcj) : 
where Jac is the Jacobian between the variables / and Therefore, 
3ac~ 1 e- H dvol = e~T, 1 ' JJdl' d(p, 

i.e. I's and ^>'s are independent, the I's are exponential and the </>'s are uniform. 
Such a decomposition was proved by McKean-Vaninsky 113. It can be useful in 
computing different thermodynamical quantities, as explained below. 

3. Thermodynamics of Moser-Calogero Potentials. 

Consider particles on the line interacting with the potential U = 2x~ 2 (the rational 
case) or 2sinh~ 2 a; (the trigonometric case). To construct the Gibbs' state for 
an infinite system of such particles, we periodize the potential as in U Ul (») = 
U{* + kbj\). In our case, U u is simply the Weierstrass 2p with possibly infinite 

k 

second imaginary period u>2 (the rational case). Now N particles on the circle of 
the perimeter ui\ are governed by the Hamiltonian 

N 2 N 
k—1 k^j;k,j—l 

The partition function is defined as 

Z{N^l3) = i| e-P H ^d N qd N p, 

M 
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where the integral is taken over the phase space M — ([0, cji) x i? 1 ) N and j3 = T^ 1 
is the reciprocal temperature. In the thermodynamic limit N and io\ tend to oo , 
but their ratio u>\/N tends to some constant: the specific volume v > 0. It is well 
knownu that, in the limit, the specific free energy per particle 

tp(v,/3) = lim N- 1 log Z(N, uj x , 0) 

exists in the limiting ensemble. The average energy per particle and the pressure 
can be computed by simple differentiation of ip 

dp ' dv 

Assuming the possibility of decompositioning for Gibbs' states described in section 
2, we compute ip(v,/3) by some version of the method of stationary phase. 

Obvious rotational symmetry allows one to reduce the dimension of the phase 
space. Indeed, according to the classical prescription, Arnold lij, 

h = -J pdq = -J^ Pk dqk = -^ Pk = -P, 

where P is the total momentum. The corresponding angle ipi is just the usual on 
the circle. 

The change of variables q k — ► q' k and pk — > p' k + jf makes the total momentum 
P' vanish. The construction of the action-angle variables Lz~tp2, • • • , In, Pn on the 
reduced phase space is based on the KP-1 equation, see E-3. We do not need an 
explicit representation for them now. 

To compute the partition function, the domain of integration is similarly re- 
duced: 

Z{N,w lt p) = fa J e-WwO / S(P(q,p) - P )dP d N qd N p 

M Rl 

= f dP faJ e^ H ^S(P(q,p) - P )d N qd N p. 

R 1 M 

The canonical transformation q^ — > q' k and p k — > p' k + t& produces 
H(q,p) = H(q',p') + ^, P = P' + P a 

and 



Z = J dP ^ J e~^'^P P o' 2N 8{P')d N q' d N p' 
R 1 M 

= V^NT^ J e^ H ^5{P) d N q d N p. 
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Now make a canonical transformation to the action-angle variables Ii,<pi,- ■ ■ , In , ip^. 
Then 



1 

m 
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Integrating out I± and all the angles we obtain 
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To proceed farther we have to make few assumptions. The arguments in favor of 
them are supplied at the end of this section. 

Assumption 1. The energy H can be expressed as H = ^2^=2 I'k ("trace for- 
mula" ) , where I' 2 , • ■ • , I' N are the actions relative to some higher symplectic structure 
to'. 

Assumption 2. The range of the variables I' 2 , ■ ■ ■ , I' N is the rectangular domain 
{(iLr ■ ■ , I' N ) ■ ^fcmin^fc' fc = 2, ••• ,7V}, as shown in the picture 



JY 



J' 

N mill 



1 2 min 



Then i7 m i n = Y] J^„ and in the thermodynamic limit, H m i n /N — > h(v), where 
h( v ) = Sfc^o U(kv) is the specific energy of one particle in the ground state with 
the specific volume v of the limiting infinite system. 

Assumption 3. The Jacobian between the variables I%, ■ ■ ■ , In and 1%, • • • , I' N 
in the thermodynamic limit has the asymptotics 

iV- 1 log Jac(4 • • • , I' N , ut, N) = logiV - s(v, 0) + o(l) 

in probability relative to the Gibbs' ensemble. The nontrivial part of the asymp- 
totics s(v,f3), we call entropy, by analogy with statistical mechanics, Lanford E3. 
Now we employ the trace formula: 

-, , N 

Z(N,cj 1 ,(3) = V2^NTcj 1 (27r) N ~ 1 — / e -"E*i Jac(7^, • • • , I' N ) J] dl' h 

J I n 

range of i' 2 ,—,i' N 

For the specific free energy in the thermodynamic limit, we obtain finally 
i/>(v,j3) = limJV _1 log Z(N,wx,P) 



k=2 
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log 2tt + N- 1 log -L / e -' 3 E^e wlo « N -^'«nd4 + o(l) 

range of i' 2 ,— ,i' 



k=2 

N 



log27r-7V- 1 logiV! + iV- 1 log J e^^^fldl'k 



range of l' 2 ,---,i' N 



+ logN_ s(v,P) +o(l) 



e 



-0H„ 



= log 2^ + 1 + A' 1 log " - s(v, /3) + o(l) 

= log27r + l-^(u)-log / 9-a(u, / 9). 
The average energy is 

<2?>=-*^=T + fc(t,) + 4(^/3)- 
The equation of state has the form 

— dv — = ~~ v ^ ~ 

These formulas provide a bridge between two unsolved problems: the asymptotics 
of the Jacobian determines the thermodynamic functions and vice versa. We will 
try to guess an explicit form of s(v,/3) in the next section. 
Now we explain why we expect our assumptions to be true. 



Assumption 1. The trace formula appears: in the 



, _ inear problem, 0, as the 

Parsevall identity and, in the cubic Schrodinger case El , as a result of contour 
integration on the Ricmann surface. We expect such formula here too, with I's 
being a periods of some Abelian differential. 

Assumption 2. A general fact due to Atiyah-Guillemin-Sternberg states that, 
for a compact symplectic phase space with a torus action, the image of the momen- 
tum map is a convex polytope. The shape of the momentum map here is suggested 
by the previous work on PDE and by the fact that the energy in the ground state 
is strictly positive. 

Assumption 3. To explain why the Jac^i ■ • • > I'n, wi, N) becomes a constant 
in thermodynamic limit, we have to define all objects under consideration on one 
probability space. Consider a new phase space M for the infinite-particle system 
on the entire line. M contains all possible positions and velocities of particles, with 
the sole restriction that any finite region in the configuration space contains a finite 
number of particles (so-called locally finite). 

Let M Ul C M comprise all spatially wi-periodic configurations. On M Wl , let 
A*wi(*) b e a Gibbs' state with periodic boundary conditions. Namely, consider 
projection of any configuration from M UJl into a finite volume, from to wi. Any 
configuration is determined by the projection. Define on the projections a Gibbs' 
state with periodic boundary condition, ft induces a measure /i Wl (») on M u>1 and 
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on the big space M, too. A limiting Gibbs' state /j. o( < ) on M is obtained from 
by passing to the limit u>\ — > oo. 

On the big space M acts the one-parameter group of spatial translations Tb,b G 
R . Obviously, and M Ul itself are invariant under the action of r. The 

limiting measure Hoo(») is also invariant under such translations and even ergodic. 

For any such M u)1 we can define a Jacobian Jac(l2, • • • , I' N , wj., N); it can be 
considered on the whole M as a function defined almost everywhere with respect to 
/ji Ul (•). The Jacobians are r-invariant functions. If they have a limit defined almost 
everywhere with respect /i o( , ) J then it is r-invariant and must be a constant due 
to the ergodicity of this measure. The exact form of the asymptotics is suggested 
by @. 

4. Alternative Approach. The Rational Case. 

Gallavotti-Marchioro and also Francoise 111, considered the system of N particles 
on the line with the Hamiltonian 

N 2 \2 2 N 
k—1 k^j;k,j—l 

They obtained a closed expression for the partition function: 

Z(N, A, (3) = ^ J e-? H d N q d N p=± (^) exp(-/3AA(7V - 1)). 

R N xR N 

Now, pass to the limit with A — > oo and assume A ~ A _1 c(?;, (3), where c(v,(3) 
is some unknown function. With this scaling iJ grow like A in the statistical 
ensemble. In the limit, we obtain an infinite-particle system with specific volume 
v and temperature (3. For the specific free energy we have 

0) = N-i log ± N e^- 1 ) + o(l) 

= - A -1 log A! + log 2vr + log A - log (3 - log c - /3c + o(l) 
= log 2-7T + 1 — /3c — log /3 — log c. 

The expression for the specific free energy contains the unknown function c(v,/3), 
but the role of the two unknown functions c(v, (3) and s(v, (3) is different. Comparing 
these equations we conjecture c(v,f3) = h(v) and s(v,(3) — logh(v). In the rational 
case h(v) can be computed explicitly 

k^O k=iO 

Therefore 

d , v , 2T 4tt 2 , , 2tt 2 

p = -T— log/it;)-^ = — + <^>=T + ^( i ;)=T+— . 
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These equations are similar to the ideal gas for which p = T/v and < E >= T/2. 
In the trigonometric case (U (x) = 2 sinh~ 2 x) we expect similar results. 
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